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1.  INTRODUCTION 
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This  paper  investigates  the  heat  budgets  produced 
by  a  limited  area  heat  source,  as  modelled  in  a  large-eddy 
simulation.  The  physical  prototype  is  a  series  of  Arctic 
leads,  breaks  in  the  Arctic  pack  ice  which  allow  relatively 
warm  ocean  water  to  contact  the  cold  atmosphere.  This 
idealized  case  represents  periodic  inhomogeneous  surface 
conditions.  Local  budgets  of  temperature  and  vertical 
turbulent  heal  flux  will  be  compared  to  their  domain 
averaged  values. 


stratification  established  upwind  of  the  lead  is  eroded 
by  the  turbulence  created  by  the  warm  surface,  producing 
a  relatively  well-mixed  region.  Below  this  region  a  stable 
internal  boundary  layer  grows  downwind  of  the  lead. 
The  figure  does  not  cover  the  entire  horizontal  domain, 
so  differences  between  the  left  and  right  borders  represent 
the  increasing  stratification  which  results  further 
downwind. 


2.  LARGE-EDDY  SIMULATION 

A  previous  large-eddy  simulation  which  modelled 
a  single  lead,  by  not  allowing  a  parcel  to  cycle  through 
the  domain  (Giendening  and  Burk.  1992).  has  been 
extended  for  approximately  three  parcel  cycles  with 
periodic  lateral  boundary  conditions,  representing 
periodic  surface  forcing  as  might  occur  over  multiple 
leads. 

Idealized  initial  conditions  consist  of  a  air-water 
temperature  difference  of  25K  over  a  200m  wide  lead, 
below  a  linear  stratification  of  56-52  =  lOK/km.  The 
large-scale  eeostrophic  wind  speed  is  2.5m/s  perpendicular 
to  the  lead  edge.  These  conditions  generate  a  large 
surface  temfjerature  flux  over  the  lead,  averaging 
0.162Km/s.  Over  the  surrounding  ice  the  surface 
temperature  flux  is  downward,  due  to  the  warmer  air 
above  it.  averaging  -0.0076Km/s.  The  heat  source  covers 
only  9%  of  the  domain,  so  the  domain  averaged 
temperature  flux  is  -t-0.0071Km/s. 

The  large-eddy  model  used  is  fundamentally  that 
described  in  Moeng  (1984),  except  that  surface  inhomoge¬ 
neity  is  introduced,  a  radiative  boundary  condition  is 
employed  at  the  model  top,  the  mixing  length  is  required 
to  meet  similarity  restrictions  near  the  surface,  and  the 
surface  similarity  includes  stable  conditions.  Further 
details  are  given  in  Giendening  and  Burk  (1992).  This 
simulation  requires  much  smaller  grid  spacing  than  that 
of  a  convective  boundary-layer  (BL).  Vertical  resolution 
is  4m  and  horizontal  resolution  is  8m.  The  grid  covers 
2.J(Mm  perpendicular  to  the  lead.  1 92m  transversely,  and 
120m  vertically  (x,y, 2  directions  respectively).  The  time 
step  is  0.5s.  Local  averaging,  both  in  time  and  parallel 
to  the  lead,  is  required  to  obtain  adequate  statistics;  the 
results  presented  are  for  t=0.6-0.7hr.  so  local  averages 
are  based  upon  17,280  values. 

3.  LOCAL  AVERAGES 

Fig.  I  gives  contours  of  potential  temperature  6, 
where  the  overbar  represents  a  local  average.  The  stable 
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F/g.  /:  Potential  temperature  6  (in  "Cj. 

Because  much  of  the  domain  is  stably  stratified,  the 
’large-eddy"  results  apply  only  to  the  region  of  the  grid 
within  which  turbulence  is  resolved.  Fig.  2.  giving 
contours  of  w  '  -  vs  here  '  represents  a  deviation  from 
the  local  average,  indicates  that  region.  For  a  parcel, 
travel  time  over  the  lead  is  relatively  short  —compared 
to  the  Brunt-Vaisalla  period-  so  the  maximum  vertical 
velocity  of  an  individual  thermal  occurs  considerably 
downwind  of  the  lead.  Here  I  use  ’thermal"  to  indicate 
an  individual  eddy,  whose  combined  effects  produce  the 
"plume"  apparent  in  Fig.  2.  Note  that  the  downwind 
tilt  of  the  plume  signifies  that  an  individual  thermal  which 
reaches  its  maximum  vertical  velocity  later/sooner  than 
the  "typical"  thermal  is  usually  both  further/closer 
downwind  and  further/closer  from  the  ground;  this  tilt 
does  not  characterize  the  individual  thermals,  which  are 
essentially  vertically  oriented. 

Fig.  3  displays  the  total  vertical  turbulent  temperature 
flux,  summing  the  resolved  and  unresolved  contributions. 
The  upward  flux  is  strongest  at  the  surface,  with  a 
maximum  axis  tilting  downwind.  The  heat  plume 
downwind  of  the  lead  results  from  the  net  heat  carried 
Ly  individual  thermals.  Downward  heat  flux  occurs  in 
the  downwind  portion  of  the  plume  and  also  at  the 
surface. 


of  the  plume  region,  but  its  magnitude  is  smaller  than 
expected  for  a  homogeneous  convective  BL.  Only  a 
minor  portion  of  the  vertical  heat  transport  occurs  through 
mean  motion,  a  consequence  of  mean  vertical  velocities 
being  relatively  small.  Note  that  the  vertical  turbulent 
flux  based  upon  perturbations  from  the  domain  mean 
is  simply  the  sum  of  the  local  mean  and  turbulent  fluxes. 
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Fig.  3:  Total  vertical  turbulent  temperature  flux 
(in  xICf^Km/s). 

Although  not  shown,  the  profile  of  w"%  where  " 
represents  a  deviation  from  the  domain  average  T~).  is 
both  qualitatively  and  quantitatively  similar  to  that 
expected  for  a  convective  BL  over  a  homogeneous 
surface,  when  non-dimensionalized  by  w.s[gQ,D/0J‘  ’ 
where  Q.  is  the  domain  averaged  surface  temperature 
flux  and  D  is  the  plume  depth  as  defined  by  the  level 
of  maximum  domain-averaged  downward  turbulent  heat 
flux.  Non-dimensionalization  of  this  inhomogeneous 
case  using  homogeneous  convective  BL  scaling  based 
upon  domain  means  is  thus  appropriate,  despite  the  large 
horizontal  inhomogeneities  and  the  limited  region  of 
strong  turbulence  downwind  of  the  unstable  surface. 
This  non-dimensionalization.  which  will  be  called 
"periodic"  scaling  to  emphasize  its  employment  under 
inhomogeneous  conditions,  allows  comparison  to  results 
obtained  for  a  horizontally  homogeneous  BL. 


4.  TEMPERATURE  BUDGETS 

Fig.  4  presents  dimensional  profiles  of  domain 
averaged  vertical  temperature  flux.  The  total  flux  profile 
is  qualitatively  similar  to  that  over  a  homogeneous  heated 
surface,  but  here  the  linear  decrease  with  height  does 
not  indicate  a  well-mixed  region,  since  such  is  not  found 
in  the  domain  mean  temperature  profile,  but  simply  that 
the  shape  of  the  temperature  profile  remains  constant 
with  time.  As  would  occur  for  homogeneous  surface 
heating,  there  is  a  net  downward  heat  flux  near  the  top 


Fig.  5:  Non-dimensional domainaverage of  b  d /dt budget 
terms:  mean  vertical  flux  -d  w  6  /dz  (dashes),  turbulent 
vertical  flux  -dw  ‘  9  '/dz  (thin  solid),  subgrid  (dots),  and 
total  (thick  solid). 

Fig.  6  presents  local  profiles  at  three  locations,  spaced 
equally  about  the  position  of  maximum  plume  growth; 
upwind  (x=400m),  at  (x=600m),  and  downwind 
( X  =  800m)  .All  perturbations  are  based  upon  deviations 
from  the  local  average,  not  from  the  domain  mean 
average. 

Upwind  of  the  position  of  maximum  plume  growth 
(Fig. 6a),  vertical  turbulent  flux  creates  cooling  at  the 
surface,  which  is  countered  by  horizontal  transport  of 
warm  air  immediately  above  the  lead.  Turbulent  heating 
occurs  above  this,  through  a  depth  about  equal  to  that 
of  the  lower  region,  an.'^  -s  countered  by  a  mean  hcruontal 
flux  which  is  now  relatively  cool.  The  remaining  terms 
are  relatively  small  in  magnitude.  The  net  3O/0t  has 
a  warming  maximum  at  0.4D. 


Fig  6'  Non-dimensional  dO/Bt  budget  terms  at:  (a) 
x=400m  (b)  x=600m  (c)  x=800m^  Symbols  as  for 
Fig.  5  plus  mean  horizontal  flux  -d  u  0  /Bx  (thick  dash) 
and  turbulent  horizontal  flux  -Bu  '  6  IBx  (dash-dot). 

At  the  position  of  maximum  plume  growth  (Fig. 6b), 
results  are  qualitatively  similar  to  those  upwind  but  are 
smaller  quantitatively  and,  in  addition,  cooling  resulting 
from  mean  vertical  transport  is  opposed  by  heating  from 
horizontal  mean  transport,  resulting  from  the  upward 
tilt  of  the  isentropes. 

Downwind  of  the  position  of  maximum  plume  growth 
(Fig.6c),  the  terms  are  further  reduced  in  magnitude. 
The  primary  qualitative  difference  from  the  previous 
locations  is  the  additional  cooling  just  above  D  resulting 
from  downward  turbulent  transport,  which  helps 
the  warming  horizontal  transport.  The  net  warming  B  6  /dt 
now  occurs  at  a  higher  level  that  at  upwind  positions. 

Locally,  therefore,  mean  horizontal  transport  is  of 
major  importance.  The  sign  of  the  horizontal  transport 
changes  with  height  -warming  near  the  surface,  cooling 


above,  and  then  warming  near  the  plume  top-  and 
opposes  the  effect  of  vertical  turbulent  transport.  The 
net  horizontal  turbulent  transport  is  negligible  at  all 
locations,  although  it  can  be  significant  for  an  individual 
eddy.  For  the  domain  average  of  the  local  terms  (Fig.5), 
the  horizontal  transport  terms  must  sum  to  zero  due  to 
the  lateral  periodicity. 

5.  VERTICALTURBULENTHEATFLUXBUDGETS 

Fig.  7  presents  domain  averages  of  the  terms  of  the 
local  Bw'  6  '  /dt  budget,  non-dimensionalized  by  the 
periodic  scaling  w.Q,/D,  for  comparison  to  the  local 
budgets  presented  below.  These  profiles  are  not 
equivalent  to  the  corresponding  terms  from  domain 
budgets  of  dw‘'d'7dt,  but  are  generally  similar  in  form 
since  mean  vertical  motion  is  small  and  local  averages 
usually  not  too  different  from  the  domain  averages.  The 
primary  difference  of  the  dwV(di  budget  lies  in  its 
buoyant  and  pressure  terms,  which  become  increasingly 
larger  near  the  surface. 


Fig.  7:  Non-dimensional  domain  average  ofBwfJ  'JBt 
budget  terms:  vertical  temperature  gradient -w 'm'  'B6/dz 
(thick  solid),  vertical  turbulent  transport  -d  '  w  '  d  '  /Bx 
(solid),  buoyant  s$  '6  '  /9„  (thick  dashes),  pressure 
correlation  -p  ' 6  Bp  '  /Bz  (dashes),  and  subgrid  (dots). 

The  profiles  in  Fig.  7  are  similar,  both  qualitatively 
and  quantitatively,  to  those  obtained  over  a  homogeneous 
surface  except  that  the  gradient  production  term  is 
significantly  larger  above  the  mid-plume  level.  Buoyant 
production  is  positive  with  maxima  at  the  ground  and 
near  D.  Its  maximum  at  D  is  of  smaller  magnitude  than 
that  of  a  homogeneous  BL.  Buoyant  production  is 
balanced  primarily  by  pressure  correlation  only  in  the 
lower  plume,  gradient  production  being  more  significant 
in  the  upper  plume.  Turbulent  transport  exports  heat 
flux  from  the  lower  to  upper  plume.  Vertically  integrated 
gradient  production  is  negative,  due  to  the  stratification 
into  which  the  plume  must  grow. 

Fig.  8  presents  local  dw  '  d  '  /dt  budgets  at  the  three 
locations  used  for  the  temperature  budgets,  for  comparison 
with  the  domain  averages  of  Fig.  7.  Several  terms  of 
a  horizontally  inhomogeneous  heat  budget  have  not  been 
plotted  because  their  contribution  is  relatively  small. 

Upwind  of  the  position  of  maximum  plume  growth 
(Fig.  8a).  buoyant  production  is  the  largest  positive  term, 
countered  primarily  by  vertical  turbulent  transport. 
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Gradient  production  is  a  significant  positive  terms,  not 
the  negative  term  of  the  domain  average.  Pressure 
correlation  counters  positive  turbulent  transport  in  the 
upper  portion  of  the  turbulent  region. 


Fig.  8.-  Non-dimensional  '  d  '  /dt  budget  terms  at: 
(a)  x=400m  (b)  x=600m  (c)  x=800m.  Symbols  as 
for_Fi^.  7_plus  horizontal  mean  flux 
-d  u  w  9  /dx  (dash-dot). 

At  the  position  of  maximum  plume  growth  (Fig. 8b), 
the  buoyant  production  maximum  has  moved  upward 
and  weakened,  as  have  the  maxima  of  positive  and 
negative  transport.  Gradient  production  is  now  negative 
throughout  the  plume  depth,  reaching  a  maximum  just 
below  D.  At  most  levels  the  negative  pressure  correlation 
makes  only  a  secondary  contribution. 

Downwind  of  the  position  of  maximum  plume  growth 
(Fig.  8c),  buoyant  production  is  still  the  primary  positive 
term  but  its  maximum  lies  Just  above  D.  The  primary 
negative  term  is  gradient  production. 


Overall,  buoyant  production  is  the  primary  positive 
budget  term,  its  maximum  moving  upward  and  decreasing 
in  magnitude  downwind.  Gradient  production  is  initially 
positive,  but  quickly  becomes  a  major  negative  term  due 
to  the  stable  stratification  which  the  atmosphere  seeks 
to  establish.  Unlike  the  temperature  budget,  the  horizontal 
transport  contribution  to  the  w  '  0  '  budget  is  small  for 
the  local  budgets.  The  local  terms  sum  to  give  the  domain 
average  effect,  which  is  much  smaller  than  local  term 
values.  Turbulent  transport  is  strongest  where  w  '  •  is 
largest,  in  the  initial  growth  region  of  the  plume.  Nec,‘ 
the  plume  top,  gradient  production  —not  pressure 
correlation—  is  the  primary  negative  term.  The  region 
downwind  of  the  plume  makes  the  largest  contribution 
to  downward  w'  6  '  . 

6.  CONCLUSIONS 

The  heat  flux  generated  over  a  limited-area  heat 
source  is  transported  vertically  by  individual  thermals 
to  create  a  mean  plume  downwind.  Domain  averaged 
terms  of  dd/dt  and  dw  '  9  '  /dt  (and  dw”d’73t)  are 
generally  similar  to  those  over  a  homogeneous  surface, 
qualitatively  and  quantitatively.  when_norvdimensionalized 
by  periodic  scaling.  For  the  dw  '  3  '  /dt  equation, 
however,  gradient  production  is  significantly  larger  than 
over  a  homogeneous  surface,  due  to  stable  stratification 
surrounding  the  plume,  and  dominates  the  pres^re 
correlation  term  near  the  plume  top.  For  the  local  d  9  /dt 
equation,  mean  horizont^  transport  is  of  major  impor¬ 
tance. 

The  similarity  of  the  domain  averaged  budget  terms 
to  those  of  a  homogeneous  convective  BL  depends  upon 
mean  transport  being  small;  for  simulations  with  larger 
mean  vertical  transport,  as  when  the  flow  is  more  parallel 
to  the  lead  and  advection  perpendicular  to  the  plume 
decreases  and  departures  occur.  Thus  the  results 
presented  here  apply  only  when  mean  vertical  motion 
is  small. 

There  are  many  caveats  for  any  simulation  such  as 
this.  1  believe  that  its  most  significant  limitations  are: 
(1)  the  results  depend  upon  the  chosen  parameterization 
for  the  subgrid  length  scale,  which  especially  affects  the 
transition  from  resolved  to  unresolved  turbulence  at  the 
plume  edge,  and  (2)  the  horizontal  resolution  is  twice 
the  vertical  resolution. 
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